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A  COMPARISON  OF  EXACT  AND  HEURISTIC  ROUTINES 
FOR  LOT  SIZE  DETERMINATION  IN  MULTI-STAGE  ASSEMBLY  SYSTEMS 

Introduction 

The  classical  or  Wilson  economic  lot  size  model  determines  the  lot  size  that 
minimizes  the  sum  of  production  and  inventory  carrying  costs  for  a  single  stage 
system.   Demand  is  assumed  to  be  continuous  and  constant  with  no  stockout  per- 
mitted.  A  number  of  extensions  to  the  basic  single  stage-model  have  been  devised 
[4  I  including  provisions  for  non-instantaneous  production,  and  for  discrete,  but 
constant,  demands.   A  further  large  class  of  extensions  considers  stochastic  de- 
mands.  The  distinguishing  characteristic  of  all  of  these  models  is  that  the  ob- 
jective is  minimization  of  costs  over  an  infinite  horizon. 

A  different  fundamental  approach  to  determination  of  lot  sizes  is  based  on 
the  assumption  of  discrete  known  demands  occurring  through  a  finite  horizon.   Such 
an  approach  allows  consideration  of  non-constant  demands  and  a  time  varying  objective 
function.   Manne J6 1 ,  Dzielinski  and  Gomoryj 3  ,  H.  Wagner  and  WhitinI  9 j,  and 
H.  Wagnerj  8  j  develop  results  for  a  single  facility.   Dantzigl2  j introduces  the 
notion  of  multi-facility  systems  where  production  of  items  at  one  facility  requires 
inputs  from  other  facilities  (with  a  linear  cost  structure) .   Veinott  10  I  and 
Zangwill  [  11,  12,  13,  14  [consider  extensions  to  concave  cost  objectives  including 
the  important  case  of  a  production  set  up  charge  with  linear  production  and  hold- 
ing costs.   The  authors  in  a  companion  paper  j  1  I  discuss  the  special  case  of  a 
finite  horizon  assembly  system,  where  each  facility  (stage)  supplies  inputs  for 
one  successor  stage. 

The  multi-stage  assembly  system  for  an  infinite  horizon  has  been  analyzed  by 
Schussel  I  7  I .   We  consider  in  this  paper  the  same  problem  with  a  slightly  modified 
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objective  functions  as  derived  in  1  1  |  .   Heuristics  for  lot  size  determination 
are  presented  and  a  dynamic  programming  algorithm  described.   Computational  re- 
sults for  all  routines  are  reported. 

The  Multi-stage  Lot  Size  Model  ' 

At  each  level  of  a  multi-stage  assembly  system  an  operation  is  performed  that 
either  acquires  raw  material,  performs  a  manufacturing  operation,  or  performs  an 
assembly  operation.   These  operations  may  take  place  at  one  geographic  location, 
for  example  in  the  intermittent  manufacture  and  assembly  of  units  in  a  job-shop 
or  flow  shop.   They  may  also  occur  in  a  multi-stage  production-warehousing  system. 
Figures  1  and  2  show  examples  of  such  multi-stage  systems. 

Let  F  =  •(  F^  ,  F„, ,  F  I  be  a  set  of  N  stages  in  a  multi-stage  system.  We 

define  the  unit  at  F  as  A  and  the  relationship  F  <'  /  F  between  stages  F  and  F 

nn  m^^n  mn 

to  denote  that  F  is  an" immediate  predecessor"of  F  .   That  is  to  say  that  one  or 
m  n 

more  units  of  A  are  used  directly  in  the  assembly  of  A  .   Similarly  F  /  F  implies 
m  ^  ^  n  ■'     m^  n   ^ 

that  F  is  a  "predecessor  of"  F  and  A   contains  one  or  more  units  of  A  . 
m       '^  n      n  m 

In  the  example  of  Figure  2,  A,  is  assembled  from  one  unit  of  A,  and  A^  and  A^ 

contains  one  unit  of  A, .   The  routines  discussed  in  this  paper  will  handle  the 

case  in  which  F  has  one  or  more  immediate  predecessors  and  contains  one  or  more 
n  '^ 

units  of  each  predecessor.   They  will  not  handle  the  case  in  which  F  precedes 
more  than  one  stage,  although  heuristics  have  been  devised  for  this  case   5 J 

At  each  stage  of  the  multi-stage  network,  F  ,  we  assume  a  fixed  set-up  cost, 
S  ,  and  a  variable  cost  per  unit,  C  ,  comprising  the  cost  of  assembly  and  incre- 
mental manufacture. 


(S>^^— <!^^^S> 


FIGURE  1 
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FIGURE  2 


Given  the  assumption  of  instantaneous  production  and  setting  Q  =  Q  , 

n=l,  2 »N-  1  ,then  the  system  contains  inventory  only  at  the  final  stage,  F  . 

N 

At  each  intermediate  stage  F  ,  Q   is  produced  and  immediately  transmitted  to 
the  succeeding  stage.   If  0  ^  Q  ,  F  <^<(  F  ,  then  in-process  inventory  occurs  at 
stage  F  .   We  define  the  relationship  between  lot  sizes  at  successive  stages  as 

k^  =  Q^  where  F^«  F^  and  K^  =  Q^   Finally  let  J^  =  juJ.  k^,  ....,  k^,^.  l} 

On  %' 

i     *       *  *   1  ) 

be  a  particular  set  of  lot  size  relationships  ,   J^  =  <  k  ,  k_ ,  ....,  k  -ijl? 

the  optimal  set, and  J  =  1^  a  set  of  all  ones. 

Note  that  we  consider  only  solutions  with  a  single  lot  size  associated  with 
each  stage.   This  rules  out  solutions  which  require  more  than  one  lot  size  at  a 
stage.   This  restriction  may  be  reasonable  in  practice  as  a  consequence  of  the 
cost  of  lot  size  changes.   Computational  experience  indicates  that  optimal  solutions 
obtained  with  this  restriction  vary  only  slightly  from  the  true  optima. 

It  is  shown  in  flj  that  the  optimal  set  of  lot  size  relations  J^  given  the  above 


* 
assumption  is  characterized  by:   (a)  k  >1 

it 

(b)  k   integer 


/    n=l,  2,  ....,N. 


These  results  are  extensions  of  those  obtained  for  the  finite  horizon  case 
8,  9,  12   .   Together  (a)  and  (b)  express  the  condition  of  optimality  that  pro- 
duction at  a  stage  occurs  only  if  production  occurs  simultaneously  at  its  successor 
stage.   Figure  3  shows  the  inventory  at  each  stage  of  a  three  stage  production 
process  with  k  =  2  and  k„  =  3,  with  the  integrality  condition  satisfied. 

It  is  possible  to  express  the  average  inventory  at  F  as  a  complex  function 

m 

of  Q  and  Q  where  F  <<F  .   The  value  of  A  is  then  taken  to  be  the  sum  of  the 
m      n        m    n  m 

variable  costs  contained  in  its  components.  The  average  value  and  the  total 

carrying  cost  can  then  be  computed  for  the  stage.   A  simpler  formulation  assumes 

each  stage  is  independent  with  the  constraint  on  positive  integer  values  for  k 

n 

imposed  in  the  search  algorithm.   Instead  of  computing  the  average  inventory  at 
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6. 

a  stage,  this  model  computes  the  average  number  of  units  in  the  complete  system 
that  have  undergone  the  activity  of  a  stage.   For  F„  of  Figure  1  we  would  calcu- 
late the  average  inventory  at  F„  and  F„  and  value  it  at  C„ ,  rather  than  calcula- 
ting the  average  at  F-  and  value  it  at  C^  +  C„. 

Since  the  demand  on  the  total  system  is  assumed  to  be  constant,  the  system 
inventory  of  the  product  of  any  stage  F  will  decline  linearly.  Given  the  assump- 
tion of  the  integrality  of  k  ,  there  will  be  exactly  k  production  lots  of  Q  , 
for  each  lot  of  Q  ,  F  <^<^F  .  At  the  point  in  time  that  Q  is  produced,  0  will 
also  be  produced  and  at  the  instant  before  that  production  begins  the  inventory 

of  A  at  F  and  A  at  F  will  be  zero.   The  argument  may  be  extended  to  all 
m     m      n     n  o  j  

successors  of  F  .   Thus  at  a  point  in  time  all  units  of  A  are  out  of  the  system 
m  m  — ^ 

at  all  stages.   At  that  point  a  lot  is  produced  and  the  system  inventory  becomes 
Q  .   Thus  the  system  inventory  of  A  has  the  familiar  saw-tooth  pattern  of  the 
simple  E.O.Q.  model.   In  Figure  4  the  system  inventory  for  the  problem  of  Figure  3 


(v^ 


is  shown.   The  cost  of  inventory  for  the  production  of  stage  F  will  be  (Q  -1)  C  I 

T 
if  discrete  demand  for  A^  is  assumed  and  where  I  is  the  carrying  cost.   If  R  is  the 

total  period  demand. the  set-up  cost  will  be  R  S   and  the  total  cost  for  A  will  be 
■  ^  ^  —  m  m 


(vvv-^^  ^ 


\ 


m  '   The  total  cost  for  the  system  give  Q^,  and  J.  is  then: 
TOTAL  COST  =   2^ 


n  =  1 


We  now  differentiate  this  function  and  solve  for  an  optimal  Q.,  •  given  J.. 


"  'n.J 


2«Efa 


Y-  then  Qjj,^   =    Fq  lorrQj+  1  where j^    Idenotes 

\^  I  y"c  K 
Z-a  n  n 

the  integer  part.     Similarly  for  the  situation  of  non- instantaneous  production 


for  P  -rP  ,  F  //F   for  F  ,  m  =  1,2,...,  N-1.   Where  P  is  the  production  rate 
mZ  n   mW  n       m  m 

K 
n 


V  C  K  /l-R  \ 


The  heuristic  solution  routine  presented  in  the  next  section  essentially  searches 

* 
on  the  elements  of  J,  then  given  a  particular  set  J.  solves  for  the  optimal  Qj^  .  and 

the  cost  associated  with  the  solution.   The  dynamic  programming  routine  is  also 

simplified  by  the  independence  of  adjacent  stages  that  has  been  presented  here. 


the  computation  of  a  lower  bound  on  the 

I  yn  "  Vc  I  +  R  s   >  >  Q„  =  .  riTT" 


This  view  of  the  problem  also  allows  the  computation  of  a  lower  bound  on  the 

N 
complete  problem  of   TOTAL  COST  >  -r— \ 

r  -J  n  =  1  T»  n 

In  reference  1  this  bound,  and  extensions  of  it  are  used  to  limit  search  for 

optimal  solutions. 

Heuristic  Search  Routines 

The  three  heuristic  rules  presented  in  this  section  search  on  the  elements 
of  k  c  J  in  an  attempt  to  reduce  the  total  cost  of  the  system.  Figure  5  shows 
the  basic  search  routine  used  in  all  three  heuristics.  To  begin  k^  ^  is  increased 


by  one,  the  K  are  updated  then  Qj,  and  total  costs  are  calculated.   If  there  is 


an  improvement  kj^,-,  is  again  incremented  and  the  process  is  repeated.   If  not, 
kjj_  is  decremented  and  a  test  of  cost  is  made.   On  the  first  pass,  given  that 
k  j^  =  1  in  the  beginning,  and  that  it  was  possibly  increased,  it  will  not  be 
economic  to  reduce  it.   However,  in  later  passes,  reduction  may  be  indicated.   The 
routine  would  then  reduce  k    until  it  was  no  longer  economic,  then  move  on  to  the 
next  element  k  _„  first  increasing  then  decreasing  it.   All  elements  in  J  are 
then  examined  and  modified  when  appropriate.   The  simplest  "single  pass"  heuristic 
stops  here.   The  second  heuristic,  "multiple  pass",  continues  to  iterate  through  J 
until  no  improvement  is  found. 


FLOW  CHART  FOR 
BASIC  SEARCH  HEURISTIC 
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The  final  heuristic,  "modified  multiple  pass",  is  identical  to  the  multiple 
pass  heuristic  except  that  it  does  not  begin  with  J  =  1.   Instead  for  each 

F  ,  n  =  1,  2,  . . . . ,  N  we  calculate  "unconstrained"  lot  sizes  Q  =  R  S 

C  I 
n 

Then  beginning  with  Q  ,  rounded  to  the  nearest  integer  value,  we  adjust  the  lot 

size  of  all  predecessors  of  F  ,  say  F  ,  to  the  integer  multiple  of  Q  nearest 

their  "unconstrained"  lot  size.   Each  predecessor  of  F  would  then  have  its  lot 

^  m 

size  adjusted  to  an  integer  multiple  of  the  adjusted  Q  and  so  on.   This  process 
defines  a  beginning  vector  J  which  may  be  closer  to  the  J^  than  a  J  =  1  would 
be.   Computational  results  for  all  heuristics  are  given  in  the  next  section. 

The  Dynamic  Programming  Algorithm 

Economic  lot  sizes  can  be  determined  by  dynamic  programming.   This  technique 

yields  the  optimal  solution  at  the  cost  of  increased  solution  time.  We  define 

the  dynamic  programming  states  to  be  the  integer  quantities  at  each  stage.   Solution 

proceeds  from  raw  material  stages  to  the  final  stage.   The  recurrence  relation  is 

defined  as  follows : 

Let«a^(n)  denote  the  set  of  predessors  of  F  for  all  m  such  that  F  //  F  . 

'^  m  m^^  n 

Let  I  denote  the  set  of  integers/  A  1.   u  (Q)  =  the  minimum  lot  size  cost  for  F 

and  all  predecessors  of  F  assuming  Q  is  the  lot  size  for  F  .  Then 

n  .^^  n 

U  (Q)  =(Q-1  )  C  I  +  R  S      2^  minimum  (  /  Q) 

"     \  2  /  "    —^     ""me^Cn)         /  _  [ 

Procedures  for  surmounting  computational  difficulties  associated  with  the  unbounded 

state  space  are  described  in  I  1  1. 


10. 

Computational  Results 

The  routines  of  the  previous  two  sections  were  tested  on  six  problems,  all 
using  the  multi-stage  network  of  Figure  2.   The  problems  are  outlined  in  detail 
in  Appendix  I.   All  programs  were  written  in  PLl  and  run  on  the  MULTICS  time- 
sharing systems  of  PROJECT  MAC,  M.I.T.,  operating  on  a  GE645. 

Table  1  illustrates  the  improved  start  that  the  "Modified  Multiple  Pass" 
heuristic  can  provide.   If  problem  two  is  solved  with  J  =  1,  the  minimum  total 
cost,  used  as  a  beginning  solution  is  $4,283.   The  improved  start  has  a  cost  of 
$2,928  and  the  heuristic  terminates  in  one  step  with  a  final  cost  of  $2,922. 
The  optimal  solution,  achieved  with  the  dynamic  programming  algorithm,  is  also 
shown. 
Sample  Run  of  Modified  Multiple  Pass  Heuristic  -  Problem  Two 


START:   J^  =  |  6, 3, 3, 4, 3,2, 4, 7, 4, 7, 3,2*2 ,2  ,3,3, l} 


Q.,  „  =  1-^    Cost  =  2,928 

N  ,0 


PASS  1:   Jj^  =1  5, 3, 3, 4, 3, 2, 4, 8, 3, 7, 3, 2, 2, 2,2,3, 1> 


Q*   =  16    Cost  =  2,922 

N,l  ' 


HALT: 


OPTIMAL  .  ^ 

SOLUTION:   J*  =  <  6,3,3,4,3,2,5,8,4, 7, 3, 2, 2, 2, 2,2,4, 1> 

Q.,  .  =  12    Cost  =  2,920 
N,* 

Table  1 
The  results  of  all  routines  on  all  problems  are  shown  in  Table  2.   Lower  bounds 
for  each  problem  are  given,  as  well  as  the  minimum  cost  with  J  =  1_.   For  these 
costs  and  network  structures  the  optimal  solution  is  close  to  the  lower  bound 
in  some  cases  and  quite  far  from  it  in  others.   The  same  is  true  for  the 
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* 
relationship  between  the  unconstrained  Q  and  Q^  ^  .   For  all  problems  the  cost 

of  the  best  heuristic  solution  is  very  close  to  the  optimal  solution  and  both 

* 
Q  ,  and  J  from  the  heuristic  are  close  to  optimal  solution.   It  is  clear  in  all 

cases  that  the  solution  derived  from  J  =  1  is  a  poor  one.   In  Figure  5  we  show 

the  relationship  between  total  cost  and  Q  for  the  case  J  =  1  and  J.. 


CONCLUSION  AND  EXTENSIONS 

Economic  lot  sizes  in  multi-stage  assembly  systems  can  be  determined  by 
dynamic  programming  for  problems  of  moderate  size,  and  by  heuristics  for  larger 
problems.   The  techniques  described  herein  can  be  generalized  to  handle  structures 
slightly  more  complicated  than  the  pure  multi-stage  assembly  case;  in  particular, 
sub-assemblies  may  be  produced  to  satisfy  exogenous  stages.   In  addition,  more 
general  concave  cost  functions  can  be  handled  provided  that  the  integrality  con- 
dition is  satisfied. 


14. 


APPENDIX  I 


PROBLEM  DATA 


15. 


H 
W 

<: 

c 

ooooooooooooooooo 
ooooooooooooooooo 

CN            i-HiHC^i-l                    CMCN                    n           iH^tH 

ooooooooooooooooo 
ooooooooooooooooo 

rHOOinr-(Ou~iO<NLnOr0cv4O"^'-(O 
rOiHiHiHcMcO                M                CS                rH 

H 
W 

<: 

Q 

OOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOO 

r-irHin.HiHino>HioojcMominooo 

r-l             1-1                     .-1                     T-l    (N    <■ 

ooooooooooooooooo 
ooooooooooooooooo 

mooiouioiAoiomoO'^O'AiriiH 

iHfO(NCMiH.H           (M           CslCMiHcsirOt-l 

r-l 

H 
W 

a 

c 

OOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOO 
OOOOOrHOOOOOOOOiHOrH 
i-lrHOiHrH           l-lr^r^l-lI-^0^-lI^           iH 

OOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOO 

rH,-IOi-lr-IOOiHOiH.HO.Hi-IOOO 
O               OOO               O               ooo 
i-t                  i-HOO                  rH                  tHOi-I 
iH           iH                                                >H 

H 

iHcsiro^invor>.oooNO'HcMmsru-ivDr^ 

iHiHiHrHiHiHrHr-t 

o 

in 

o 

CVJ 

o 

• 

f-l 

II 

II 

pi 

l-l 

16. 


H 
W 

<; 

H 

<: 


w 
<: 

H 

<: 


H 

w 

< 

< 
p 


OOOOOOOOOOOOOOOOO 

ooooooooooooooooo 


C>J    t-t    -*    in    O     rH 


o  o  o 

CN    rH    CSI 


o  o  in  csi 


u 


ooooooooooooooooo 
ooooooooooooooooo 

ininiOLominmmmLomLniommmio 


u 


ooooooooooooooooo 
ooooooooooooooooo 

LTiiTimiriminiriLoinmLniniriu-iLrimu^ 


ooooooooooooooooo 
ooooooooooooooooo 

mi-<mosiom>— ifHtooocNOiooioin 

^  iHrH  i-lrH.H<NiHi-l 


ooooooooooooooooo 
ooooooooooooooooo 

OOOOOOOmomiriu-iinOLriOio 


OOOOOOOOOOOOOOOOO 
OinoOOOOOOOOOOOOOO 


r-irHLorocMroomOfMrH^d-rni-l 


in  lO  in 


w 

H 

0:1 


rHcNro-vt-invDr^oooNOr- tr-jro<rinvor~ 


o       m 

O  O^ 

o 


»J  M 


17, 


REFERENCES 


1.  Crowston,  W.  B.  and  M.  Wagner,  "Multi-Stage  Assembly  Systems,"  M.I.T. 
Working  Paper  (forthcoming). 

2.  Dantzig,  G.  B.,  "Optimal  Solution  to  a  Dynamic  Leontief  Model  with 
Substitution,"  Econometrica  23,  1955,  pp.  295-302. 

3.  Dzielinski,  B.,  and  R.  Gomory,  "Optimal  Programming  of  Lot  Sizes,  In- 
ventories, and  Labor  Allocations,"  Management  Science,  Vol.  11,  No.  9 
(July  1965),  pp.  874-890. 

4.  Hadley,  G.  and  T.  M.  Whitin,  Analysis  of  Inventory  Systems,  Prentice- 
Hall,  Inc.,  Englewood  Cliffs,  N.  J.,  1963. 

5.  Henshaw,  A.  C. ,  "Multi-Stage  Inventory  Models ,"  M.I.T.  S.M.  Thesis, 
1969. 

6.  Manne,  A.  S.,  "Programming  of  Economic  Lot  Sizes,"  Management  Science, 
Vol.  4,  No.  2   (January  1958),  pp.  115-135. 

7.  Schussel,  G, ,  "Job-Shop  Lot  Release  Sizes,"  Management  Science  Vol  14 
No.  8  (1968),  pp.  B449-B472.  ^ '        ' 

8.  Wagner,  H.  M. ,  "A  Postscript  to  'Dynamic  Problems  of  the  Firm,'  " 
Naval  Res.  Logistics  Quarterly,  7,  pp.  7-12  (1960). 

9.  Wagner,  H.  M. ,  and  T.  M.  Whitin,  "Dynamic  Version  of  the  Economic  Lot 
Size  Model,"  Management  Science,  Vol.  5,  No.  1,  October  1959,  pp.  89-96. 

10.  Veinott,  A.  F. ,  Jr.,  "Minimum  Concave-Cost  Solution  of  Leontief  Sub- 
stitution Models  of  Multi-Facility  Inventory  Systems,"  17,  Operations 
Research,  March  -  April  1969,  pp.  262-291. 

11.  Zangwill  ,  W.  I.,  "Models  of  a  Dynamic  Lot  Size  Production  System," 
Management  Science,  Vol.  15,  No.  2,  May  1969,  pp.  506-527. 

12.  Zangwill  ,  W.  I.,  "Minimum  Concave  Cost  Flows  In  Certain  Networks," 
Management  Science,  14,  1968,  pp.  429-450. 

13.  Zangwill  ,  W.  I.,  "A  Deterministic  Multiproduct  Multifacility  Production 
and  Inventory  Model,  Operations  Research ,  14,  1966,  pp.  486-507. 

14.  Zangwill.,  W.  I.,  "A  Deterministic  Multi-Period  Production  Scheduling 
Model  with  Backlogging , "  Management  Science,  13,  1966.  pp.  105-119. 


m 


t^-fw 


££C2    >r 


USFM^ 


^.;t 


Date  Due 


P.?- 


*^^ 


^i^ 


iHkf 


?  i  Jic 


wsa 


APR.  2  3  1935 


jaw: 


f^.\0^ 


8Cti  ^  T 


1^ 


19  » 


Lib-26-67 


t//-'--  7,7 


3  ^060  003  TOb  HDO 


Toa 


■Kr-^ 


D  DQ3  fi75  Sbl 


3  '^DfiO  003  675  5T5 


3  T060  003  TO 


E  555 


3  TDaO  003  TOb  531 


MIT  LIBRARIES 


3   TOfiO   0 


03 


b    bl4 


•^70-7  0 


Hyh  70 


MIT  LIBRARIES 


3  TOfl 


0  0 


03  fl75  E 


^/:^-7<:^ 


3  TOflO  DD3  fl7S  SE 


Q71-70 


MIT  LIBRARIES 


</7'^-  70 


3    TDflO    003    TOb    5MT 


W7 


J- 70 


TOaO    003    675    b03 


